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ABSTRACT 



Context. Within low-mass star formation, water vapor plays a key role in the chemistry and energy balance of the circumstellar ma- 
terial. The Herschel Space Observatory will open up the possibility to observe water lines originating from a wide range of excitation 
energies. 

Aims. Our aim is to simulate the emission of rotational water lines from envelopes characteristic of embedded low-mass protostars. 
A large number of parameters that influence the water line emission are explored: luminosity, density, density slope and water abun- 
dances. 

Methods. Both dust and water emission are modelled using full radiative transfer in spherical symmetry. The temperature profile is 
calculated for a given density profile. The FLO level populations and emission profiles are in turn computed with a non-LTE line code. 
The results are analyzed to determine the diagnostic value of different lines, and are compared with existing observations. 
Results. Lines can be categorized in: (i) optically thick lines, including ground-state lines, mostly sensitive to the cold outer part; (ii) 
highly excited (E u > 200 - 250 K) optically thin lines sensitive to the abundance in the hot inner part; and (iii) lines which vary from 
optically thick to thin depending on the abundances. Dust influences the emission of water significantly by becoming optically thick 
at the higher frequencies, and by pumping optically thin lines. 

Conclusions. A good physical model of a source, including a correct treatment of dust, is a prerequisite to infer the water abun- 
dance structure and possible jumps at the evaporation temperature from observations. The inner warm (T >100 K) envelope can be 
probed by highly-excited lines, while a combination of excited and spectrally resolved ground state lines probes the outer envelope. 
Observations of H' s O lines, although weak, provide even stronger constraints on abundances. 
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1. Introduction 

Low-mass (M <3 M ) young stellar objects (YSOs) form 
through gravitational collapse of cloud cores. Most of their mass 
is accreted over pe riods of less than 10 6 years dAndre et al.[ 2000; 
Mvers et al. 2000). These early stages of low-mass star forma- 
tion, the so-called Class and Class I stages, are character- 
ized by the presence of centrally con densed envelopes which 
veil t he central protostars and disks dLadalfl98 7; Andr e et ail 
119931) . The envelopes also contain most of the mass (M ~0.1- 
1 M Q ) of the total system until the emergence of the central star 
and d isk dAdams et al.ll9 87; Shirl ey et al.l20 00; J0rge nsen et al.l 
|2002|) . Water plays a pivotal role in these early evolutionary 
stages. First, the high dipole moment of H2O makes its ro- 
tational lines excellent coolants of warm gas in the inner re- 
gion of protostellar envelopes, all owing material to continue 
collapsing at higher temperatures dGo ldsmith & Langerl 119781: 
iNeufeld & Kaufmari[l993t ICeccarelli et al.lll996l) . Second, wa- 
ter is one of the dominant r eservoirs of the non-refra ctory oxy- 
gen budget in dense clouds dvan Dishoeck et al 1 fl993l The par- 
titioning of gaseous oxygen between H2O, OH and O, as well 
as the ratio between H2O gas and ice, affect the chemistry of all 
other o xygen-containing mo lecules including the complex or- 
ganics (Charnle v"eTaIlfl992h . 



The evapo ration point of H 2O ice in interstellar space is 
around 100 K C Fraser et al.ll200ll) . so only the small warm inner 
regions of protostellar envelopes are likely to contain fractional 
abundances of water vapor as high as 1 0" 4 with respect to H2 , 
perhaps even exceeding that of CO (e.g.JCeccarelli et al1 ll996; 
Ivan D ishoeck & Blake 1998; Boonman et al.l2003l) . In colder re- 
gions, water is the main ingredi ent of icy grain mantles as shown 



by infrared spectroscopy (e.g.. iWh ittet et al. 



1989; Boogert et al. 2002; Pontoppidan et al 



19831; ISmifh et al.l 
2004), with orders 
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of magnitudes lower gas-phase water abundances. This abun- 
dance jump makes water an excellent tracer of warm, dense gas 
and an unique probe into the inner region close to the forming 
star. 

Observations of interstellar gas-phase water from the ground 
are limited by the large amounts of water vapor in the Earth's 
atmosphere. Water maser lines at radio wavelengths long-ward 
of 3 mm (100 GHz), such as at 22 GHz , are readily seen , 
even from low-ma ss embedded YSOs (e.g.. lFuruva et al.| [2003; 
iPark & Choil 120071) . but are constrained to small hot-spots and 
often variable in nature dFelli et al.ll2007l) . They provide little 
information about the general physical an d chemical evolution 
during star formation dElitzur et al. 1 119891) . One of the few ex- 
ceptions of a non-masing water line observable from the ground 
is the 183 GHz J = 3o - 2o2 transition, which can be used to 
trace wide-spread water emission under excep tional conditions 
(e.g. 



wiae-spreaa water emission under excep' 
. ICernicharo etanfT99Ql [19961: lDot^l2000h . 
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Observations of deuterated water can also be done from 
the ground, in particular through the HDO lot - O po 464 GHz, 
3 12 — 2 21 225 GHz and 2n — 2i2 241 GHz lines (e.g., Turner et al. 
1971 iJacaet all Il990t iHelmichet al] 11991 IStark et al.1 12004: 
Parise et al l 120051) and the D 2 ho - loi 316 GHz line 
dButner et al.l 120071) . but the large uncertainties in the D/H ra- 
tio owing to fractionation make such observations a poor tracer 
of water itself. The H^O isotopologue can be detected as 
well, but only through its high excitation tem perature (E up = 
193 K) 3 n -2?o transition at 203 GHz (e.g.. iJaca etalJ [1988; 

1.1 120061) . 



1998 
200C 



iGensheimer et al.ll 19961; Ivan der Tak etldBooi 

To observe water emission lines systematically for a wide 
range of energy states, one has to turn to space-based ob- 
servatories. Gaseous water has been detected in a number 
of star forming regions using the Infrared Space Observatory 
(ISO), the Sub-millimeter Wave Astrono my Satellite ( SWAS) 
and ODIN (e.g. iHelmich et all 119961: lHarwit et al.l 
Nisini etalJ 1 19991: ICeccarelli et al.1 Il999t ISnelletai] _ 
Berein & Sneill 120021: iMaret et al.l 120021: INisini etaD l2002t 
Boonman etaLll2003l: iBergin & Melnickll2005l: iRistorcelli et all 
2005; ICernicharo & Crovisierl 120051) . However, due to large 
beams, low spectral resolution or a combination of both, such 
observations were unable to accurately resolve the water abun- 
dance profile and trace the origin of the water emission, es- 
pecially for low-mass YSOs. Two instruments on-board the 
Herschel Space ObservatorjQ can be used for future observations 
of rotational water line s at unprecedented se nsitivity, resolution 
and spectral coverage (ICeccarelli et al. 200 1|). The Heterodyne 
Instru ment for the Far Infrared (HIFI) ( de Graauw & He lmich 
|2001|) a nd Photodetector A rray Camera and Spectrometer 
(PACS) dPoglitsch et al.ll2001l) provide high spectral (up to 125 
kHz) and spatial (20" at 1 THz) resolution, which will allow 
a thorough analysis of water emission through a much wider 
range of observable transitions. In addition, the high sensitiv- 
ity of the instrument will allow observations of various lines at 
low abundance, such as the optically thin isotopologues of wa- 
ter, in particular H^O. However, to infer abundances from such 
observations, detailed radiative transfer calculations are needed. 

Several modelling efforts have been carried out to analyze 
the observed water lines, either b ased on models of shocks 
at high and low velocities (e.g ., Kaufman & Neufeld 1996; 
iGiannini et al.ll200U INisini et al 112002b or on thermal emission 



from the envelope itself (e.g., Maret et al. 2002; Boonman et al. 
I2003S) . As pointed out most recently bv lvan der Tak et al.l (12006). 
the high water abundance combined with its high dipole moment 
and its susceptibility to far-infrared emission from dust results in 
high optical depths and make this molecule particularly difficult 
to model. As a result, no systematic parameter studies have yet 
been carried out. We present here a parameter study of the water 
emission from low-mass protostellar envelopes using a full ra- 
diative transfer method including the effects of dust and applied 
to realistic physical models of low-mass YSOs. Specific predic- 
tions for Herschel are made by convol ving the computed emis- 
sion w ith the Herschel observing beams. Poelma n& van der Takl 
d2007l) present similar models for high-mass protostars. 

This paper is organized as follows. §2 presents the physi- 
cal and molecular parameters necessary to accurately model wa- 
ter. §3 takes an in-depth look at a single Class source, L 483, 
and explores the effects of possible water abundance profiles. 
§4 presents the results from the range of parameters that influ- 
ence water emission. §5 investigates the influence of dust on wa- 
ter line fluxes and profiles, which is often ignored. §6 gives an 



analysis of a selected sample of observed water lines from ISO, 
SWAS and ODIN to test the models, and to infer implications 
for observations with future instruments such as HIFI and PACS. 
Final conclusions and recommendations can be found in §7. 

2. Model 

The modelling has been carried out in two stages. First, the 
dust radiative transfer is calculated for a given power-law den- 
sity profile (§2.1). This produces dust temperature profiles for 
given density structures. Subsequently water level populations 
and line emission profiles are mo delled using the line radia- 
tive t ransfer code RATRAN (§2.2) dHogerheijde & van der Takl 
2000), with the assumption that the gas temperat ure is equal to 
that o f the dust. For an overview of the steps, see Ivan Dishoeckl 
(2003) (Figure 2). Our method refers to the 'Empirical method' 
presented there. See TableQ]for an overview of the parameters. 

2.1. Physical structure of the envelope model 
2.1 .1 . Dust radiative transfer 

Our protostellar envelope models are taken to be spherically 
symmetric, with the protostar located at the center. It has been 
shown that this assumption is valid for the analysis of both 
the continuum as we ll as line emission down to scales of a 
.IShirlev et alJl2OOOl;|j0rgensen et al.|[2 002; 



few hundred AU (e . 
ISchoier et 



jd AU (e .g 

IDISoqIIj 



J0rg ensen et al.ll2005l) . The density n at ra- 



dius r within the protostellar envelope is described by: 
n(r) = n (r/r y p , 



(1) 



see http://herschel.esac.esa.int/ 



with the power law exponent p and density no at a reference ra- 
dius ro as parameters. Within the context of protostellar collapse 
theories, both parameters depend on the age of the system and 
on the original core conditions, such as its initial mass. 

Both the gas and dust temperatures inside the envelope are 
often also described by a power-law profile. However, this ap- 
proximation fails at small radii where the dust becomes op- 
tically thick (e.g.. | Dotv & Leungl 1 19941 : |j0rgensen et al.l [2002; 
YoungjfcEvanl2005). 

The actual temperature structure depends 
on a combination of the density structure and the central source 
luminosity, and a dust radiative transfer calculation is needed to 
determine it. The radiation is produced by the energy released 
through accretion of material onto the protostar. 

The continuum radiative transfer through the envelope was 
calculated using th e spherically symmetri c dust radiative trans- 
fer code DUSTY dlvezic & Elitzurf [19971) . DUSTY is a scale- 
free 1-D code, and uses the ratio Y of the outer over inner ra- 
dius, R ut/R\n, the dust opacity t at 100 fim, the temperature at 
R m and the power law index of the density profile p as its vari- 
ables. Results are then scaled by using the luminosity and dis- 
tance of specific sources. The inner edge of our envelope R m was 
set at the radius where the temperature of the envelope reaches 
250 K. Material closer to the star is assumed to be located in 
a disk, which is not taken into account. In addition, we assume 
that the disk has no influence on the temperature and density 
distribution of the envelope. The exact inner radius in AU de- 
pends on the luminosity, which differs for each model. R^ is 
typically of the order of a few tens of AU, much smaller than 
the Herschel beams (see Table[2]). In all our models Y is taken to 
be 200. Values of ^ out are a few thousand to ten thousand AU, 
larger than the Herschel beam (20" at 1113 GHz) at all but the 
lo west frequencies , but g enerally smaller than the values given 
in |j0rgensen et al.l d2002l) . The adopted dust opacity's are the 
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often-used OH tables, column 5 of Ossenko pf & Henningl (fl994) 
throughout. Figure Q] shows a typical temperature and density 
profile. 

2.1.2. Gas radiative transfer 

The gas temperature is taken to be equal to the dust tempera- 
ture for the entire spatial extent of our model, and the gas:dust 
mass ratio is taken to be 100:1. Warm dense gas can efficiently 
cool through the emission of ex cited molecular and a tomic lines 
such as CO, H 2 0, OH and [O I] dGiannini et al]200ll) . Water can 
also heat the gas through fa r-infrared pumping fo llowed by col- 
lisional de-excitation ("e.g.. ICeccarelli et alJI 19961) . As a result, 
dust and gas temperat ures can differ slightly, especially in th e 
100-200 K range (e.g JCeccarelli et al.ll2000tlMaret et alj|2002l) . 
However, at densities of 10 7 cirT 3 or higher the gas and dust in 
the warm inner region couple efficiently. At low densities (<10 5 
cm ~ 3 ) and temperatures (<50 K) any differences between gas 
and dust temperature depend on the source parameters (e.g., lu- 
minosity, total e nvelope mass) and the ex ternal interstellar radia- 
tion field (ISRF). lDotv & Neufeldld 1997b show in their Figure 12 
that these differences are generally small, of order 10 K or less 
for low- mass YSOs, and typical ly only exist at radii larger than 
2000 AU. lBoonman et al J 020031) show that for high-mass YSOs, 
the water line emission for a full radiative transfer calculation of 
the gas temperature differs negligible from that calculated using 
the 7dust = r gas assumption. 

Our choice of /? out =200 Ri n assumes that water vapor at 
larger radii does not contribute to the total water line emission 
due to the low densities (< 10 5 cirT 3 ) and temperatures (<15 
K) in these regions. This assumption mostly affects the transi- 
tions leading to the ground states. Any such material is able to 
absorb water emission for optically thick lines, and lead to weak 
emission lines in observations taken with a beam that is much 
larger than our adopted source sizes. However, such extended 
envelope material is almost indistinguishable from any larger 
scale cloud complexes not directly related to the source. Such 
contributions from large-scale clouds need to be taken into ac- 
count on a source-by-source basis, but are beyond the scope of 
this work. The protostellar envelopes are assumed to be shielded 
from strong external radiation fields by surrounding molecular 
cloud material. If directly exposed to outside radiation, the tem- 
perature in the outer envelopes is affected by the ISRF, but ma- 
terial is not heated by more th an 40 K for the ISRF f ound within 
the Orion star-forming region dj0rgensen et al.l2006l) . Low-mass 
star-forming regions are expected to have significantly lower 
ISRF. External radiation is not considered for these models. 



Table 1. Parameters for low-mass protostellar envelopes 



Variables 



;iH 2 (1000AU)(10 b cm- j ) 


0.4, 1, 5 


P 


1.5,2.0 


Luminosity (L G ) 


2, 7, 25 


X (water) 


10~ 4 , 10~ 6 


Xi (water) 


10" 6 , 10" 7 , 10" 8 


Fixed Parameters 


Distance 


150 pc 


Stellar temperature 


5000 K 


Gas temperature 


= Tdust 


R m (inner radius) 


250 K 


R out (outer radius) 


200 R m 


ro 


1000 AU 


Water freeze-out 


100 K 


o : p ratio F^O" 


3:1 


o : p ratio H 2 


1:1 


Dust properties 


OH 5 b 


Velocity dispersion AV 


1.67 km s-' 



"For the L483 models, ortho- and para-water were modelled indepen- 
dently. See §3 6 See Ossenkopf & Henning (1994) column 5 

Table 2. The masses and inner radii of the model envelopes. 





«H 2 (ro) 








L e 


10 6 cm" 3 


M G 


AU 


AU 


p = 1.5 


2 


0.4 


0.06 


8.22 


23 


2 


1 


0.15 


8.22 


23 


2 


5 


0.76 


8.22 


23 


7 


0.4 


0.15 


15.3 


42 


7 


1 


0.39 


15.3 


42 


7 


5 


1.93 


15.3 


42 


25 


0.4 


0.40 


29 


80 


25 


1 


1.01 


29 


80 


25 


5 


5.04 


29 


80 


p = 2 


2 


0.4 


0.07 


8.22 


21 


2 


1 


0.18 


8.22 


21 


2 


5 


0.90 


8.22 


21 


7 


0.4 


0.13 


15.3 


38 


7 


1 


0.33 


15.3 


38 


7 


5 


1.65 


15.3 


38 


25 


0.4 


0.24 


29 


72 


25 


1 


0.63 


29 


72 


25 


5 


3.15 


29 


72 



2.1.3. Grid parameters 

To limit the model complexity, the large range of parameters was 
reduced to four variables (see TableQ}. The density at 1000 AU, 
«o(H2), and the exponent of th e density power law p are two 
of these. J0rge nsen et alj d2002l) observed many low-mass YSOs 
and inferred densities at 1000 AU between 1 x 10 5 to 5 x 10 6 
cirT 3 . Accordingly, three values, 0.4 x 10 6 , 1 x 10 6 and 5 x 10 6 
cm -3 were chosen to sample this range. Similarly, p was chosen 
to be either 1.5 or 2.0. The stellar variables are reduced to a sin- 
gle variable, the luminosity. Models using luminosities of 2, 7 
and 25 L are explored to represent the large range of luminosi- 
ties observed in low-mass protostars. Although the surface tem- 
perature of the protostar cannot be derived directly, it is assumed 
to be around 5000 K for low-mass young stars, regardless of the 



luminosity or age. Test models show that due to the reprocessing 
of stellar photons into IR/thermal photons very close to the star, 
the stellar effective temperature has little influence on the mod- 
elling. The final variable is the gas-phase water abundance, to be 
discussed below. TableQ]summarizes all the adopted parameters. 
The resulting masses and inner and outer radii for each envelope 
model can be found in Table [2] 

2.2. H 2 line modelling 

Because of the gradient in temperature, the water abundance 
X with respect to H2 is expected to vary within YSO en- 
velopes. Previous studies of low- and high-mass YSOs used a 
variety of H 2 observations with ISO-LWS, ISO-SWS and/or 
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Fig.l. FIGURE TOO LARGE FOR ASTRO- 

PH. CAN BE DO WNLOADED FROM WEBSITE 
http://www.strw.leidenuniv.nl~kempen/water.php Example 
of the physical structure of a model envelope with the density 
(solid), temperature (dashed) and abundance (dotted) displayed 
versus radius. The adopted parameters are «o(H2)=10 6 cm" 3 , 
p=1.5 and L — 2 L . 



SWAS to constrain the water abundance profile in protostel- 
lar envelopes to have a 'jump' in abund ance at a charac- 
teristic temperature a round 100 K (e . g.. | Maret et all 120021; 
iBoonman et al] 120031: Ivan der Tak etai] 120061) . Such 'jumps' 



have also been confirmed by ground-based obs ervations of other 
molecules such as CH^OH a n d H?CO (e.g. , van der Tak et al 
2000; CeccaiellTetai] 120001; ISchoier et al] 



sumed to be in a 1:1 ortho to para ratio. Water itself is in 3:1 
ortho to para ratio. Collisional rate coefficients for both ortho- 
and para-H^O with ortho- and para-H2 were obtained from the 
L AMDA datab a se (|Schoier et al] 120051) based on calculations 
by iGreen et al.l dl993 ) for temperatures in the range from 20 
to 2000 K. New rate coefficients have recently been published 
(Phill ips et alTl 996; D ubernet & Grosieanll2002t iGrosiean etai] 
l2003l:lFaure et alj20 07fl but those new calculations often do not 
include a sufficiently large range of temperatures nor the higher 
excited energy levels, necessary for the analysis of YSOs. Test 
models using different published rate c oefficients show that ab- 
solute line fluxes may differ up to 30% ( Poelma n & van der Takl 
120071) . Even larger differences are found at low temperatures 
(T < 20 K) and densities («(H.2) < 10 4 cm" 3 ) for the groun d 



2004; Ma ret etalj |2004). The abundance of water in the inner 
region Xq is high for two reasons. First, most water exists in the 
gas phase above 100 K, the temperature at which water ice evap- 
orates from the grains. Second, gas-phase reactions of atomic 
oxygen with molecular hydrog en drive all oxygen in to water at 
high temperatures (T>230 K) dCharnlev et al] 120011) . However, 
in regions where the temperature drops below 100 K, water is 
mostly frozen out onto the grains and only trace amounts of 
gaseous water are present with an abundance X^. 

Accordingly, 'jump' abundances are assumed for all models 
with Xq, the abundance of the warm inner region, taken to be 
10~ 4 or 10" 6 , and X^, the abundance in the cold outer region, 
taken to be 10" 6 , 10" 7 or 10~ 8 . These values cover the range 
of inferred abundances from observations to date for embed- 
ded protostars. Although beam averaged H2O ab undances below 
10~ 8 have been observed for B 68 and Oph-D (Ber gin & Snelll 
|2002|) . such clouds are starless. At the extremely low temper- 
atures of <15 K virtually all heavy elements are frozen out. 
Abund ances as low as 10~" were tested by Boo nman et al] 
(2003), but models with outer abundances of 10~ 8 were preferred 
in fitting the combined SWAS and ISO-LWS data. Note that one 
combination of abundances is a constant abundance at 10~ 6 , in- 
cluded to investigate the differences between 'jump' and con- 
stant abundance profiles. A finer grid of abundances within this 
range is explored for the L 483 model. 

Figure Q] includes the assumed abundance structure. Our 
models do not cons ider a so-called 'drop' abundance profile (see 
Schoi er et al] 12004). in which the abundance in the outermost 
part is higher due to longer timescales for freeze-out at lower 
densi ties and due to photodesorption of ices (Bergin & Melnick 
120051) . Such an abundance profile would result in more self- 
absorption in the ground-state lines, equivalent to absorption 
caused by cold fore-ground material not associated with the pro- 
tostellar system. 

The effects of micro-turbulence within the envelope are rep- 
resented by a velocity dispersion FWHM, AV=1.67 km s _I 
(b = 1 km s~ 1 , with b the Doppler parameter), consistent with the 
typical line widths observed for o ptically thin rotational emis- 
sion i n low-mass YSO envelopes (J0rgen sen et al]|2002l [2004, 
120051) . For a more thorough discussion of the effects of the 
choice of this value, see § 5.2. An infall velocity profile was 
included, based o n the velocity profiles for collapsing clouds 
as found by IShul d 1 977b . This velocity profile was scaled to a 
velocity of 4 km s _1 at the inner radius R m . At radii of a few 
hundred AU, the micro-turbulent velocity becomes dominant 
over the infall velocity, whic h is the scale probed mostly by 
single-dish observations (e.g., J0rgensen et al. 2004). H2 is as- 



state line intensities and profiles Poelman & Spaans (2005); 
Dubernet et al] d2006l e.g.,). However, the general trends found 
20021 |j0rgenseni in this paper remain valid for different rate coefficients. . For 
a general discussio n of uncertainties in molecular data, see 
Schoier et al] d2005l) . 



Table 3. 

PACS. 



Key water lines observable with Herschel HIFI and 



Transition 


Freq. 


Wavelength 


F 

"Up 


Beam" 




(GHz) 


w 


(K) 


(") 


HIFI 


Ortho-H20 transition 


1 10- lot 


556.93607 


538.66147 


61.0 


39 


2l2"loi 


1669.90496 


179.65094 


114.4 


13 


221 "2[2 


1661.00802 


180.61322 


194.1 


13 


3 12-221 


1153.12682 


260.16219 


249.4 


19 


3l2"3o3 


1097.36505 


273.38213 


249.4 


21 


321 ~3l2 


1162.91187 


257.97312 


305.3 


19 


Para-H20 transition 


1 1 1 -Ooo 


1113.34306 


269.45872 


53.4 


21 


2o2 - 1 [ 1 


987.92670 


303.66625 


100.8 


21 


2i 1 -2 02 


752.03323 


398.91854 


136.9 


30 


220 -2 11 


1228.78902 


244.14281 


195.9 


19 


331 "4()4 


1893.68681 


158.42113 


410.4 


13 


4 22 -3 3 i 


916.17164 


327.44956 


454.3 


25 


4 22 -4i3 


1207.63890 


248.41863 


454.3 


19 


5 2 4-4 3 i 


970.31524 


309.17787 


598.8 


21 


PACS 


Ortho-H20 transition 


2l2"loi 


1669.90496 


179.65094 


114.4 


13 


221-1 10 


2773.97691 


108.14798 


194.1 


8 


2 2 i-2i2 


1661.00802 


180.61322 


194.1 


13 


3o3"2l2 


1716.76979 


174.74678 


196.8 


12 


4i4-3 03 


2640.47434 


113.61595 


323.5 


8 


Para-H20 transition 


220-1 11 


2968.74895 


101.05267 


195.9 


7 


3l3"2o2 


2164.13230 


138.62371 


204.7 


10 


4fl4-3l3 


2391.57290 


125.44046 


319.5 


9 



" taken from http://herschel.esac.esa.int/home.shtml 

The excitation and radiative transfer of the rotational water 
lines were calculated u sing the spherically symmetric v ersion 
of the RATRAN code dHogerheiide & van der fakll2000l) . First, 
level populations were computed for all H2O levels with energies 
up to 2000 K for the model envelope using full Monte Carlo ra- 
diative transfer. Individual iterations use a convergence criterion 



2 The rate coefficients from lFaure et al.l ( 2007) only became available 
most of the models had been run 
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of 10~ 6 on the relative error of the fractional population for the 
solution of the statistical equilibrium. Iterations continue with 
increasing amounts of photons, until the largest relative differ- 
ence between fractional populations over all cells and all levels 
of three subsequent iterations is less than 1/6. In practice this 
means that the vast majority of cells and levels are converged 
to much smaller relative errors. Far infrared radiation from the 
dust and CMB are included in the H2O excitation and line for- 
mation. This inclusion makes it possible to investigate both in- 
frared pumping, self-absorption of a line and absorption against 
the dust continuum. Even though many levels are very sparsely 
populated, tests show that it is necessary to include the higher 
levels to obtain an accurate population distribution throughout 
the model. 

The populations are then used to reconstruct detailed line 
profiles through ray-tracing, which can be convolved with 
any beam or area. It can also be used to make velocity 
maps of the source. The adopted spectral resolution is 0.05 
km s' 1 . RATRAN was bench-marked against many codes 
( Ivan Z adelhoff et al. 2002). Furthermore, it w as part of a recent 
bench mark tesj^l for line modelling of H2O (Ivan der Tak et alj 
l2005h . 

Although level populations were calculated for all levels up 
to £u P =2000 K, only lines with an upper energy level below 500 
K are selected for the presentation of results (see Table©. Both 
H2O and H^O were modelled with a focus on lines that fall into 
the Herschel wavelength range. The 16 0: 18 iso tope ratio is as- 
sume d to be 550:1, the ratio in the local ISM (Wils on & Rood! 
119941) . All generic models are placed at 150 pc, the typical dis- 
tance of most nearby low-mass star forming regions such as 
Ophiuchus, Taurus and Chameleon. Results are convolved with 
the appropriate Herschel beams and given in main beam temper- 
atures. For reference, the typical sensitivities expected with HIFI 
(5cr, 0.5 hr, 0.5 km s 1 resolution) are 60 mK for 557 GHz, 180 
mK at 1 130 GHz and 600 mK at 1669 GHz. For PACS, limiting 
line fluxes (5<x, 0.5 hr, point source) range from 6xl0~ 18 W irT 2 
to 35 xl0~ 18 W irT 2 Q. To convert from main beam temperatures 
Tmb to fluxes S v in Jy (10~ 26 W irT 2 Hz ) use the following 
formula. 

S v = 2.65 x 10 8 x r MB (K) x 2 (') x A~ 2 (jjm) (2) 

This assumes a source size equivalent to the beam size. For 
sources smaller than the beam, a correction needs to be taken 
into account. The Herschel beam 0q ranges from 40" (0.66') at 
557 GHz (539 /mi) to 7" (0.12') at 3000 GHz (100 jum) (see 
Table©. 

3. Example : L 483 

3.1. Physical structure 

As an illustrative mo del, the isolated Class source L 483, 
(Mott e & Andrei l200 id located at 200 pc, was modelled to in- 
vestigate the sensitivity of the water line intensities to differ- 
ent abundances . This sourc e has been studie d extensively by 
J0rge nsen et al.l d2002l) and lJ0rgensenl d2004l) and its physical 
and chemical structure is well constrained through continuum 
and line data from (sub)millimeter single-dish telescopes and 
interferometers. The well-determined envelope and source pa- 
rameters (Lboi = 9 L and «o(H2)=lxl0 6 cirT 3 ) are roughly 

3 see http://www.sron.rug.nl/~vdtak/H20/ 

4 Calculated with Hspot, see http://herschel.esac.esa.int/ 



in the middle of our parameter range, making L 483 an excel- 
lent test-case to model water emission. Only the inferred density 
slope, 75=0.9, is somewhat smaller than that of other sources, but 
this does not influence the model trends. For CO, a jump abun- 
dance at an evaporation temperature of T ~30 K was found to 
explain the observed molecular lines. Any compact disk is con- 
strained to a mass of 0.012 M u sing interferometer ob servations 
with the Sub-Millimeter Array dJ0reensen et al.ll2007f) . negligi- 
ble compared to the total estimated envelope mass of 4.4 M out 
to the 10 K radius of 10,000 AU. In addition, this source has 
been observed with ISO-LWS dGiannini et al.ll200ll) . The best- 
fit parameters from the dust radiative transfer (F=1400, p-0.9, 
Tioo=0.3), which fit both the SED and the spatial extent of the 
SCUBA continuum data, determine the physical structure of the 
envelope. R m (at T=250 K) is 9.9 AU for this rather shallow 
density profile and R oM is ~ 10,000 AU. The radius where wa- 
ter freezes out onto the grains (7=100 K) is located at 35 AU. 
This L 483 model was used to explore a wider variety of wa- 
ter abundances than the generic models. Specifically values of 
1 x 10~ 4 , 5 x 10~ 5 , 1 x 10~ 5 , 5 x 10~ 6 and 1 x 10~ 6 were used for 
X and 1 x 10^,5 x 10~ 7 , 1 x 10~ 7 ,5 x 10~ 8 and 1 x 10~ 8 for 
Xj. All other parameters were kept fixed. Ortho- and para-water 
were modelled separately, and water abundances are given for 
ortho- or para-water individually instead of the combined total 
water abundance with the used 3:1 ortho to para ratio. This al- 
lows one to probe the ortho to para ratio independently. Thus, Xq 
and Xd refer to either ortho- or para-H20, depending on the line 
in question. 

3.2. Results 

In Table|4]the integrated emission of both Hl 6 and H^O lines, 
observable with Herschel-HIFI (see Table |3), are presented for 
the reference grid of L 483, whereas Table [5] presents the peak 
temperatures. Detailed contour plots of all HIFI lines as func- 
tions of Xq and Xd are presented in Fig. [2] the online appendix. 
Many of the trends discussed b elow are similar to those found 
for AFGL 2591, as discussed in Poelma n & van der Tiki (2007) 
Overall, lines are found to belong to three categories. First, all 
excited H} 8 and several H2O highly excited {E up >200 K) lines 
are completely optically thin. They have Gaussian profiles with 
widths depending on the turbulent width. Due to the freeze-out 
of water onto grains below 100 K, these optically thin lines probe 
the inner warm dense region. For the model where water has a 
constant abundance of 10~ 6 , emission from the warm dense inner 
region still dominates and is optically thin (either due to the low 
population of the highly excited states or the low abundance of 
H^O). A good example is the H 2 5 2 4-4 31 970 GHz line shown 
inFig.[2](n'g/if)- 

The ground state lines of both Hl s O and H2O, together with 
most H2O excited lines with energies below 200 K, provide a 
strong contrast with the optically thin lines. These lines only 
probe the outer region. Even an abundance of water as low as 
10~ 8 in the cold outer envelope produces a high enough column 
for these lines to become optically thick. The extremely high 
optical depth (t > 25) and the absorptions into the dust fur- 
ther prevent the integrated intensity from probing the entire en- 
velope. Thus, these lines trace the outer abundance X^, and self- 
absorption and absorption of the dust continuum are commonly 
found in their line profiles. Note, however, that most optically 
thick lines have wings that are either optically thin or less opti- 
cally thick. These wings may be able to probe the inner region 
up to the boundary layer around T = 100 K. See the 1 m-lot line 
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Table 4. Integrated intensities, J Tyy&dV [K km s _1 ], for H2O and H^O lines for L 483 models with various abundances. Negative 
values indicate strong (self-) absorbed lines. 



Transition 




X = 10- 4 






X = 10- 5 






X = 10- b 




x d = 


10- 6 


10-' 


10- 8 


10- 6 


10-' 


10- 8 


10- 6 


10-' 


10- 8 


Ortho-H 2 transitions 


1 10-I01 


2.7 


2.4 


0.83 


2.8 


1.0 


8.8(-2) 


9.5 


3.8 


0.57 


2 12- lot 


-0.2 


-0.13 


0.27 


-0.35 


-0.38 


-0.33 


3.0 


0.23 


-0.05 


2 2 i-2i2 


1.7 


0.27 


0.20 


0.17 


0.32 


0.11 


1.7 


0.37 


0.12 


3 12-221 


0.8 


0.17 


0.10 


0.47 


0.12 


7.8(-2) 


0.95 


0.20 


7.1(-2) 


3 12"3o3 


0.8 


0.29 


0.13 


0.43 


0.22 


6.0(-2) 


0.94 


0.39 


74(-2) 


321"3l2 


1.3 


0.35 


0.11 


0.73 


0.26 


5.4(-2) 


1.5 


0.46 


8.2(-2) 


Para-H 2 transitions 


1 1 1 -Ooo 


-6.5(-3) 


1.23 


-0.17 


-0.12 


5.2(-2) 


0.29 


3.12 


1.4 


-0.13 


2o2" 1 11 


1.1 


0.59 


0.26 


1.4 


0.6 


0.45 


0.6 


1.4 


0.26 


2i 1 -2o2 


0.85 


0.73 


0.16 


0.99 


0.6 


0.34 


0.96 


1.1 


0.15 


220-2n 


0.38 


0.61 


0.18 


0.45 


0.50 


0.31 


0.38 


0.83 


0.13 


331"4o4 


4.4(-2) 


6.3(-2) 


6.3(-2) 


1.9(-2) 


1.9(-2) 


2.7(-2) 


3.7(-3) 


1.903) 


1.3(-3) 


422-331 


5.0(-2) 


6.8(-2) 


3.7(-2) 


5.2(-2) 


5.6(-2) 


5.1(-2) 


1.8(-2) 


1.5(-2) 


9.7(-3) 


4 22 -4i 3 


0.21 


8.6(-2) 


5.3(-2) 


0.25 


8.3(-2) 


8.0(-2) 


0.24 


0.10 


5.8(-3) 


5 2 4-4 3 i 


6.0(-2) 


8.8(-2) 


5.5(-2) 


3.2(-2) 


4.3(-2) 


5.5(-2) 


4.0(-3) 


5.9(-3) 


3.3(-3) 


Ortho-H^O transitions 


1 10-I01 


-3.3(-2) 


-5.6(-3) 


2.3(-2) 


-4.0(-2) 


-2.0(-2) 


5.2(-3) 


-1.9(-2) 


-2.4(-2) 


-3.9(-4) 


2l2"loi 


-0.13 


-3.3(-2) 


1.2(-2) 


-0.13 


-4.5(-2) 


-6.7(4) 


-0.12 


-4.9(-2) 


-5.2(-3) 


2 2 i-2i2 


4.1 (-3) 


4.2(-3) 


3.7(-3) 


1.8(-3) 


9.6(-4) 


9.3(4) 


1.0(-3) 


1.3(-4) 


1.1(4) 


3 12-221 


5.6(-2) 


6.0(-2) 


5.9(-2) 


3.6(-3) 


3.0(-3) 


3.0(-3) 


5.9(4) 


1.39(-4) 


1.2(4) 


3 12"3o3 


2.5(-2) 


2.6(-2) 


2.5(-2) 


8.5(-3) 


6.9(-3) 


6.8(-3) 


2.0(-3) 


4.6(-4) 


3.8(-4) 


321-3 12 


1.7(-2) 


1.9(-2) 


1.8(-2) 


4.5(-3) 


3.6(-3) 


3.5(-3) 


1.0(-3) 


2.K-4) 


1.7(4) 








Para-H 


' s O transitions 










In -Ooo 


-0.21 


-7.7(-2) 


8.7(-3) 


-0.16 


-8.8(-2) 


-4.7(-3) 


-0.21 


-9.3(-2) 


-1.2(-2) 


2o2 - 1 1 1 


5.3(-2) 


4.0(-2) 


3.8(-2) 


4.4(-2) 


2.0(-2) 


1.8(-2) 


2.5(-2) 


4.7(-3) 


3.1(-3) 


2i r -2o2 


2.5(-2) 


1.9(-2) 


1.9(-2) 


2.0(-2) 


7.9(-3) 


7.6(-3) 


9.6(-3) 


1.5(-3) 


l-K-3) 


220-2 11 


1.7(-2) 


1.4(-2) 


1.4(-2) 


1.6(-2) 


6.5(-3) 


6.2(-3) 


7.2(-3) 


1.0(-3) 


6.9(-4) 


33I"4o4 


9.2(-5) 


1.2(-4) 


1.204) 


5.7(-6) 


4.6(-6) 


4.6(-6) 


5.1(-7) 


2.5(-7) 


2.3(-7) 


4 22 -3 3 i 


6.3(-4) 


8.4(-4) 


8.504) 


2.8(-5) 


2.3(-5) 


2.3(-5) 


9.1(-7) 


6.7(-7) 


6.4(-7) 


4 2 2-4i 3 


2.5(-2) 


3.0(-2) 


3.002) 


1.5(-3) 


1.2(-3) 


1.2(-3) 


4.8(-5) 


3.5(-5) 


3.4(-5) 


5 2 4-4 3 i 


7.1 (-5) 


1.1 (-4) 


1.1(4) 


1.6(-6) 


1.3(-6) 


1.3(-6) 


4.6(-8) 


4.2(-8) 


4.0(-8) 



In this and subsequent tables, the notation A(B) indicates A x 10 s . 



Fig. 2. FIGURE TOO LARGE FOR A STRO-PH. CAN BE DOWNLOADED FROM WEBSITE 

http://www.strw.leidenuniv.nl~kempen/water.php. Contour plots of the integrated water line intensity for the L 483 models, 
plotted as functions of and Xq. left : the 0-H2O 1 10-I01 li ne (557 GHz) , showing a strong dependence on X<j; middle : the p-H^O 
2o2- 1 1 1 line (995 GHz), showing a clear cut between the dependency on X& or Xq, depending on the abundance profile; right : the 
P-H2O 524-431 line (970 GHz), showing a strong dependence on Xq. 



in Fig.^(left), Tableland the spectra in Fig.|3]and|4](discussed 
in §3.4). 

Finally, some lines are optically thick for a high Xq-X^ com- 
bination, but are optically thin at low abundances. It is not known 
a-priori if these lines originate within the cold outer envelope or 
the warm dense inner regions. Examples are several of the highly 
excited H2O lines, but such dependencies are also seen for lowly 
excited H^O lines. Fig.^2\(middle) shows such an example with 

the p-H^O 2o2-l 11 line. It can be seen that, depending on the 
total column of water, this line either traces X^ at high columns, 
or Xq at lower columns. For the H^O lines a clear dividing line 
between dependency on Xq or X4 lies at X ( ;(H' 8 0) < lCT 9 5 (see 
online appendix). For H^O such a division is less clear. Note 
also that the dust can become optically thick, especially at high 
frequencies, so the higher frequency lines will not be able to 
probe the inner region at all. For the L 483 models, this only 



takes place at frequencies higher than 1500 GHz, but see §5 for 
further discussion. In the following, individual H2O and H^O 
lines are discussed in more detail. 



3.3. Integrated emission 
3.3.1. H 2 

The main isotopologue of water, Hj 6 0, is optically thick in most 
transitions considered here. All lines with J up = 1 and 2 are opti- 
cally thick and trace the outer cold envelope. Most other lines 
trace either the warm inner or cold outer envelope. The inte- 
grated intensities of these 'hybrid' lines depend sensitively on 
the precise abundance parameters of the L 483 reference model. 
In particular, the radius where the line turns optically thick is 
comparable to that where the temperature reaches 100 K for this 
model. These are the 221 - 2n, 3 12 - 221, 3 12 — 3o3, 32i - 3 12 and 
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Fig. 3. H 2 spectra for an abundance X of 1(T 5 for the 1 10 -l m (557 GHz, left), 2 12 -l i (1669 GHz middle) and 3 12 -3 03 (1097 GHz 
right) lines. Lines are shown for X& of 10~ 6 (solid), 1CT 7 (dot) and 10~ 8 (dash). 
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Fig. 4. H 2 spectra for an abundance X d of 10~ 7 for the 1 10 -l i (557 GHz left), 2 12 -l i (1669 GHz middle) and 3i 2 -3 03 (1097 GHz 
right) lines. Lines are shown for Xq of 10~ 4 (solid), 10~ 5 (dot) and 10~ 6 (dash). 



4 22 - 4i3 lines. At low abundances of X&, they depend more on 
X(). If the abundance in the outer region is larger than 5xl0~ 8 , 
they depend on X d instead. The integrated intensities range from 
several hundreds of mK km s in absorption to a few K km s _1 
in emission, with peak brightness from 100 mK to 1 K, read- 
ily detectable with HIFI. The group of highly excited para-lines 
have peaks up to a few tens of mK, requiring much longer inte- 
gration times. 

Dependencies of the line strengths on abundances are not 
linear. An exception is if the line is completely optically thin 
over the entire profile. For most lines this is not the case and 
parts of the line profile will be optically thin and parts optically 
thick, causing complicated dependencies. For instance, the to- 
tal integrated intensity of the 557 GHz line depends sensitively 
on the inner and outer abundances. Abundance combinations of 
10~ 4 /10~ 6 produce lines with less total emission than 10~ 6 /10~ 6 . 
The 10~ 6 /10 -6 combination contains less water emission origi- 
nating from the warm inner regions, but more dust continuum is 
able to escape the inner regions in the 10~ 6 /10~ 6 . This influences 
the population distribution in the outer region, resulting in less 
central absorption. Since absorption is a major feature within this 
line profile, the total integrated line intensity effectively goes up. 

3.3.2. H^O 

Because of the low abundance of the H} 8 isotopologue, most of 
its lines are optically thin. Only the three transitions connecting 
to the ground levels of ortho- and para-Hl 8 0, In — Ooo, 1 10 — loi 
and 2 12 - loi, have optical depths higher than 5 for all abun- 
dances. However, excited lines with J up —2 are optically thick in 
their line centers for Xq > 10 . The integrated intensities of 
these H^O lines range from a few hundred mK km s in ab- 
sorption to a few tens of mK km s _I in emission, whereas peak 



temperatures range from a few mK to a few tens of mK, thus 
requiring long integration times. Due to the low abundance, the 
dependence on either Xq or X& is much more evident than for 
Hj 6 (see Fig|2]for example). The optically thick ground-state 
lines accurately trace X^, while optically thin excited lines de- 
pend only on Xq. The few optically thick excited lines trace Xq 
as well, since the r = 1 surface lies within the 100 K radius. 
Higher excitation Hj 8 lines are not very strong, in the range of 
a few micro-K to a few milli-K, too weak to be detectable with 
HIFI in reasonable integration times. 

3.4. Line profiles 

In Figures [3] and @] a number of characteristic line profiles are 
presented. The continuum has been subtracted so that the profiles 
can be readily compared on the same scale. 

The total profile contains valuable information about the 
envelope structure and velocity profile, none of which can be 
determined using the integrated line intensity alone. As seen in 
the examples, there is a large variety in model profiles depend- 
ing on a number of parameters. First, the excitation energy of 
the upper state, the total water column density and the critical 
density of the line in question will determine if an emission 
line is optically thick, thin or a combination of both. Most 
lines show a combination, with optically thin wings around an 
optically thick line center. Second, the width depends on the 
systematic and turbulent velocity. Third, water in the colder 
outer envelope will absorb photons and scatter the radiation, 
causing an absorption feature. Finally, the optical depth of 
dust itself blocks line emission at frequencies above 1500 GHz 
originating in the inner region. In the following these profile 
shapes are discussed in more detail. 
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Table 5. Peak temperatures [K] of H2O lines for L 483 models of various abundances. Negative values indicate primarily absorption 
into the dust continuum. 



Transition 




X = 10- 4 






X = 10- 5 






X = 10-" 




x rf = 


10- 6 


10-' 


10- 8 


10- 6 


10-' 


10- 8 


10- 6 


io-' 


10- 8 


Ortho-ILO transitions 


1 10-I01 


2.8 


2.1 


0.65 


2.8 


1.1 


0.2 


5.4 


2.8 


0.5 


2i2-loi 


0.45 


0.35 


0.25 


0.3 


-0.13 


-0.13 


1.9 


-0.13 


0.25 


2 2 i-2i2 


1.1 


0.1 


0.05 


0.2 


0.1 


3.8(-2) 


1.6 


0.21 


4.5(-2) 


3 12-221 


0.28 


6.5(-2) 


2.4(-2) 


0.18 


4.5C-2) 


2.1(-2) 


0.36 


0.1 


2.4(-2) 


3]2-3o3 


0.27 


0.12 


4.1 (-2) 


0.12 


0.1 


1.5(-2) 


0.4 


0.17 


2.8(-2) 


321"3l2 


0.5 


0.18 


4.K-2) 


0.33 


0.13 


1.5(-2) 


0.55 


0.26 


3.2(-2) 


Para-H20 transitions 


1 1 1 -Ooo 


1.5 


1.3 


-9(-2) 


0.45 


0.45 


0.45 


1.6 


1.25 


0.13 


2o2"lll 


1.6 


0.55 


8(-2) 


1.1 


0.4 


0.18 


0.9 


1.0 


8(-2) 


2n-2o2 


0.8 


0.22 


7(-2) 


0.38 


0.18 


0.16 


0.45 


0.38 


7.5(-2) 


220-2] 1 


0.4 


0.28 


7.5(-2) 


0.12 


0.24 


0.17 


0.1 


0.35 


7(-2) 


331"4o4 


2.8(-2) 


2.2(-2) 


2.2(-2) 


l(-2) 


l(-2) 


1.9(-2) 


2.5(-3) 


1.5(-3) 


8(-4) 


422-331 


2.8(-2) 


2(-2) 


l(-2) 


2.3(-2) 


2.4(-2) 


1.5(-2) 


1.2(-2) 


l(-2) 


6(-3) 


4 22 -4i 3 


0.16 


2.4(-2) 


1.4(-2) 


0.12 


2.7(-2) 


1.7(-2) 


0.12 


4.5(-2) 


2.2(-2) 


5 2 4-4 3 i 


2.2(-2) 


2.6(-2) 


1.6(-2) 


1.8(-2) 


1.9(-2) 


2.4(-2) 


2.5(-3) 


3.5(-3) 


2(-3) 



The profiles of optically thin lines are Gaussian and are good 
tracers of the velocity profile. Optically thin lines always trace 
the inner warm region, where the higher excitation states are 
populated. The ratio of the line wing over the peak tempera- 
ture can vary, depending on the excitation energy of the line. 
Higher excitation lines are broader, while low excitation lines 
have a narrower peak with very weak line wings. An example in 
Figure|4]is the 3i2-3o3 1097 GHz line, which is completely opti- 
cally thin forXn=l()- 6 . The line center becomes optically thick at 
Xo=l()- 5 and thus decreases, but the optically thin line wings are 
noticeably higher in emission. At Xo=10~ 4 , the optically thick 
line center is brighter than at 10 due to the high abundance, 
but still weaker than the optically thin emission at Xo=10~ 6 . 

The profiles of optically thick lines are much more compli- 
cated. Their wings can either be optically thick or thin, depend- 
ing on the abundances and frequency of the line. This gives rise 
to three different profiles. 

- Lines with a Gaussian line center, but with enhanced line 
wings. A characteristic 'bump' can be seen at velocities 
where the line profile changes from optically thick to thin 
when going to more extreme velocities. The line center is 
optically thick, but is heavily pumped by the line radiation 
coming from water in the warmer regions further inwards in 
the envelope. This behaviour is often called 'effectively op- 
tically thin' . Such lines trace the inner abundance Xq, unless 
optically thick dust veils the inner region. An example is seen 
by the 3 I2 - 3 03 1097 GHz line with an X of 10" 4 in Figure 
S 

- 'Flat top' lines, where the center of the line is optically thick, 
but not self-absorbed. The center of the line is dominated by 
emission from a small region with a specific temperature. 
Line wings are optically thin. Material at larger radii is un- 
able to cause self-absorption and there is not enough material 
deeper within the envelope to make such a line 'effectively 
optically thin' . Optically thick dust can prevent warmer re- 
gions from influencing the line at higher frequencies. An ex- 
ample is given by the 3i2 — 3o3 1097 GHz line with Xj = 10~ 6 
in Figure [3] 

- Deeply (self)-absorbed lines with most of the emission in the 
line wings. Absorption into the dust occurs, and is stronger 
for high frequency lines, where the dust continuum is higher. 



All three ground state lines show such profiles. Examples are 
the lin-lni and 2i2-loi lines in Fig.[3]and[4] Both the depth 
and the width of the absorption vary, as do the line wing 
strengths: 

(i) Absorption depth: If dust is optically thick (e.g., at 1669 
GHz), the depth is independent of both Xo and X& as long 
as t >> 1, and depends only on the radius where the dust 
turns optically thick, setting the continuum level. In Figure 
[3] and [4] (middle) the 1669 GHz continuum-subtracted line 
is always at -0.2 K. If dust is optically thin (e.g., the 557 
GHz 1 10-1 01 line), the depth v of the absorption depends on 
both Xj and Xo and is due to self-absorption of the water 
in the cold outer region against the warm emission. Indeed, 
it is found that water at around 70-90 K acts as a warm 
background emitter and water at 20 K as the cold absorber. 

(ii) Absorption width: again there is a difference between 
high and low optical depths of dust. For 1669 GHz the 
width of the absorption depends on Xj, but not on Xq. 
For 557 GHz, there is no dependence on abundance. The 
width of the absorption is the same for varying Xo (the 
temperature at which the water emits changes little), and 
even though the profiles appear different for different Xd, this 
can be attributed to the changing absorption depth, not width. 

(iii) Line wing strength: the strength of the line wing, which 
often contains the bulk of the line emission, is determined by 
the optical depth of the dust and gas. Most often, the peaks 
seen in the wings of the ground-state lines are still optically 
thick. The wings trace both Xd (see 1 10-I01 an d 2i2-loi lines 
in Fig. [3] left, middle), which still contain large columns of 
the gas, as well as Xo (lin-loi and 2i2-loi Fig. [4] left, mid- 
dle), although in the latter case the relation is not straight 
forward. Line pumping from the warm gas further inwards 
is the main origin of the dependency from Xq=\0~ 4 to IO -5 . 
However, if Xd and Xo are low enough the line wings turn 
optically thin and are significantly enhanced (Fig [3] dashed 
lines) because higher temperatures are probed. Such be- 
haviour is seen for a XJXq combination of 10 _7 /10 6 where 
the line wings are brighter than for XJXq combinations of 
10- 7 /10- 4 and IQ^/IQ- 5 . 
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3.5. HIFI diagnostic lines 

Although the beams for HIFI and the accompanying dilution fac- 
tors are large compared to the physical size of the warm gas (see 
Table O, the inner warm region is detectable through the higher 
excited states. The outer cold envelope is completely transparent 
for such lines. Considering the expected sensitivities of HIFI, 
a combination of different lines can constrain the water abun- 
dance profile. The H 2 s O 2 02 -lii 994 GHz line in HIFI Band 4 
is mostly optically thin and traces Xo, except at H^O X^ abun- 
dances higher than 5xl(T 6 . However, the line is weak (a few 
tens of mK) and needs at least a few hours of integration to 
be detected. A possible alternative is the high excitation H^O 
3i2-3o3 line in HIFI Band 4, which can be detected in roughly 
an hour of integration. At high Xo abundances it traces Xo, but 
for Xq lower than 10~ 5 , emission coming from the colder outer 
part dominates, thus tracing X^. This emission originates in the 
region with temperatures ranging from 70-100 K. Once the in- 
ner abundance Xq is well-constrained, ground-state lines can in 
turn be used to constrain the abundance in the outer region. The 
bright wings of the 557 GHz line (HIFI Band 1) are detectable 
within a few minutes, but the entire line profile is needed for 
proper analysis (see §3.4). The dependence on X^ can be traced 
most accurately using the line profiles of the ln-Ooo 1113 GHz 
para-H20 and 1 101 GHz para-H^O lines in Band 4, but resolv- 
ing these profiles requires much longer integration times than the 
557 GHz line. The 1667 GHz 2i 2 -l i ortho-H 2 line in Band 6 is 
intrinsically brighter, but will take almost an hour of integration 
due to the low sensitivity of HIFI at this frequency. In addition, 
dust can become optically thick at this high frequency so that the 
line profile only traces the structure in the outer envelope. 

The optically thick H* 6 2 02 - 1 i i 988 GHz line, with its opti- 
cally thin H^O counterpart are the best candidates to determine 
the optical depth of water, but the H^O lines requires signifi- 
cant integration time (see above). An alternative is provided by 
comparison between of the 548 GHz H' s O and 557 GHz H' s O 
ground-state lines, but Hj, 8 is still optically thick so that de- 
tailed modeling is required. The ln-Ooo 1113 GHz, 2o2-ln 988 
GHz and 2n-2o 2 752 GHz para-H 2 lines are well-suited to 
probe the physical structure of the envelope: due to their dif- 
ferent excitation energies and optical depths, these lines probe 
different environments within the protostellar envelope and their 
line ratios are not affect by the uncertain ortho/para ratio. 

3.6. Comparison with ISO-LWS data 

L 483 was observed with ISO-LWS dGiannini et al.ll200ll) . The 
para 2 2 o-ln line at 2968 GHz (101 ^m) was detected with a 
flux of 8.6+1.0 x 10" 20 W cirT 2 and the ortho 2 12 -loi line at 
1670 GHz (179.5 //m) was found to have an upper limit of 4.3 X 
10 -20 W cirT 2 . The low LWS resolving power of A/AA = 200 (~ 
1500 km s _1 ) was insufficient to determine the effect of any self- 
absorption, or to disentangle outflow components from quiescent 
envelope material. The best fitting LVG model for the high-/ 
CO lines gave a component of gas with 7=850-1800 K; it is 
expected that such hot gas co ntributes significantly to the water 
emission (Gia nnini et al.|[200~Tb . 

The above models, corrected for the ISO beam of 80", are 
compared to the observed intensity. Our models show a strong 
absorption into the dust for both lines with optically thin line 
wings. Integrated intensities range from -7 x 10~ 22 up to 17 x 
10- 20 Wcm- 2 forthe2 20 -ln line and from 0.9 x lO" 20 W cirT 2 
to 16 x 10~ 20 W cirT 2 for the ortho 2i 2 -l i line for different X Q - 



X& combinations. The best fit, using both the detection and upper 
limit, is found for an outer abundance of 10~ 7 and an inner abun- 
dance between 10~ 6 and 10~ 5 . The detected line is reproduced 
within 30%, while the upper limit is only 50% higher than the 
intensity produced by the model. For lower values of of order 
10~ 8 , the integrated emission in our model is almost an order of 
magnitude lower than the observed 2 2 o-ln line. Such low outer 
abundances are f ound for cold regions in the envelopes around 
high-mass YSOs (Boonman et al. 2003), but appear inconsistent 
with the data for this low-mass YSO, unless the detected emis- 
sion originates fully in then outflow. Higher inner abundances 
cannot be ruled out. However, the self-absorption in the models 
becomes stronger, resulting in a worse fit to the detected line and 
upper limit. 

Other water lines within the ISO-LWS domain were not de- 
tected. The derived upper limits were found to be at least a fac- 
tor 5 higher than fluxes predicted from all possible models. At 
short (> 80 fim) wavelengths, the difference is a few orders of 
magnitude. Since all lines are unresolved, no conclusion can be 
reached whether the detected water emission is dominated by 
outflow contributions or by thermal excitation within the circum- 
stellar envelope with a high abundance. This example illustrates 
the need for spectrally and spatially resolved data of a number 
of lines to disentangle the different physical components. 

4. General parameter study 

Tableland Figures |5]and|6]present the results for a grid probing 
the parameters from Table Q] for a constant abundance combi- 
nation of 10~ 4 for Xq and 10 7 for X^. Results for the full grid, 
including all probed abundances, can be found in the online ap- 
pendix. The main variables of the grid - the luminosity, the den- 
sity and the density profile steepness - are discussed within this 
section. The dependence on abundance mirrors the behaviour 
seen for the L 483 model discussed in §3. For most combina- 
tions of envelope parameters, the H^O transitions studied here 
are optically thick, sometimes with optically thin line wings. For 
Hj 8 lines, most transitions are optically thin, although the three 
ground state lines are optically thick at line center. The only op- 
tically thin Hi 6 O lines are those for which the upper level of 
the transition has an excitation temperature above 200-250 K, 
specifically the lines with J — 3 and higher. These lines are 
only optically thin if the densities are low enough, <10 7 cm" 3 , 
and abundances are low, Xq < 10" 5 and < 10~ 7 . For typical 
low-mass envelopes, such densities are only encountered within 
the inner region, where enhanced abundances are encountered. 
Thus, because many higher excited H^O lines become optically 
thick only in the very dense warm inner region, they are still 
good tracers of the inner region (see §3.3.1 and Fig.[4]right). For 
high densities and thus higher envelope masses, the larger dust 
column plays a major role by heavily influencing the water line 
emission. This is discussed in §5. 

4.1. Luminosity 

Fig. [5] shows four lines for various luminosities. For all lines, 
both the line strength as well as the absorption features are 
stronger for higher luminosities. For optically thin lines, a best 
fit can be found with a slope that has a power law index of ~0.8, 
i.e. / oc L 8 . Fig.|7]shows two such power laws for the 2 2 o-2n 
1228 GHz line that differentiate between the exact conditions 
and correspond to 0.45L 8 and 0.15L 8 . The slope of the power 
law is determined by the temperature profile, which determines 
the amount of warm gas present in the dense inner region, and 
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Table 6. Integrated intensities, Jt^b^V [K km s l ], for H2O lines in the wide parameter grid for an abundance combination of 
Xq = 10" 4 and X& = 10" 7 . Column 2 to 4 show the luminosity variation with a density of 10 6 cm" 3 . Column 5 to 7 show the variation 
of density with a luminosity of 7 L 



Transition Integrated Intensity ( f T MB A V [K km s ']) 

Lum. (w. no = 10 b cm~ J ) Dens. (w. L = 7 L Q ) 

2 7 25 0.4 1 5 

p= 1.5 



Ortho-H 2 transitions 



lio-loi 


0.37 


0.95 


0.70 


0.43 


0.95 


2.1 


2l2"loi 


0.35 


0.36 


-1.2 


0.30 


0.36 


-4.4 


2 2 i-2 12 


0.29 


0.87 


1.5 


0.60 


0.87 


0.74 


3l2"3o3 


0.21 


0.60 


1.5 


0.38 


0.60 


0.89 


3 12-221 


0.21 


0.67 


1.6 


0.40 


0.67 


0.84 


3 2 i-3i2 


0.37 


1.1 


2.7 


0.59 


1.1 


1.7 


Para-H20 transitions 


I11-O00 


0.13 


1.0 


0.31 


0.43 


1.0 


1.6 


2o2"lll 


0.58 


1.5 


3.0 


0.72 


1.5 


1.1 


2n-2o2 


0.34 


1.0 


2.0 


0.41 


1.0 


2.7 


2 2 o-2n 


0.31 


1.0 


2.3 


0.47 


1.0 


1.4 


331"4o4 


4.3(-2) 


0.15 


0.58 


0.18 


0.15 


1.3(-3) 


422-331 


4.K-2) 


0.13 


0.48 


0.13 


0.13 


9.6(-2) 


4 2 2-4l3 


9.5(-2) 


0.26 


0.74 


0.18 


0.26 


2.3(-3) 


5 2 4-4 3 i 


5.0(-2) 


0.15 


0.57 


0.16 


0.15 


8.7(-2) 








P 


= 2.0 






Ortho-H 2 transitions 


1 10-I01 


0.12 


0.24 


0.19 


0.10 


0.24 


0.82 


2l2"loi 


0.89 


-1.6 


-3.1 


0.71 


-1.6 


-5.0 


221 "2l2 


0.13 


0.14 


0.15 


0.13 


0.14 


-1.4 


3l2"3o3 


0.12 


0.47 


1.3 


0.35 


0.47 


0.51 


3 12-221 


0.13 


0.48 


1.2 


0.29 


0.48 


0.49 


321"3l2 


0.20 


0.87 


2.1 


0.45 


0.87 


1.2 


Para-H20 transitions 


I11-O00 


0.30 


9.4(-2) 


0.23 


0.52 


9.4(-2) 


9.0(-2) 


2o2"lll 


0.65 


1.4 


3.1 


0.90 


1.4 


1.8 


2n-2o2 


0.46 


0.94 


2.1 


0.56 


0.94 


2.2 


2 2 o-2n 


0.32 


0.78 


2.1 


0.54 


0.78 


1.3 


331 "4o4 


-3.7(-3) 


-6.7(-3) 


0.13 


7.K-2) 


-6.7(-3) 


-1.2(-3) 


4 22 -33i 


1.0(-2) 


5.9(-2) 


0.31 


9.3(-2) 


5.9(-2) 


2.3(-3) 


4 2 2-4i 3 


4.7(-2) 


0.19 


0.74 


0.22 


0.19 


8.1(-2) 


5 2 4-4 3 i 


8.3(-3) 


6.5(-2) 


0.34 


0.11 


6.5(-2) 


-6.9(-3) 




-10 -5 5 10 -10 -5 -10 -5 5 10 -10 -5 5 10 

VKL. offset [km/s] VKL offset [km/s] VUL. offset [km/s] VEL. offset [km/s] 

Fig. 5. Examples of the luminosity dependence of line profiles for a model with p=1.5, Xq=10~ 4 , Xd=l0~ & , «o(H2)=10 6 cm 4 . From 
left to right are shown the ln-Ooo 1113 GHz, 2o2- 1 1 1 988 GHz, 2n-2o2 752 GHz and 33i-4o 4 1893 GHz lines. Luminosities are 2 
(solid), 7 (dot) and 25 (dash) L . These lines reflect the different characteristic line profiles discussed in § 3.4. It can be seen that 
the line shapes stay the same with increasing luminosity, while the line strengths and depths of the absorption are stronger at higher 
luminosities. 



the filling factor of this warm gas within the beam. This filling 
factor is not readily apparent, since models have different values 
for R m and different radii where 7=100 K is reached depending 
on luminosity and density profile (see Table 2). The peak tem- 
peratures of these optically thin lines naturally follow a similar 



dependence. Even for lines which display an optically thick line 
center (e.g., the 33i-4o4 line in Fig. [5j, the peak temperatures de- 
pend on luminosity in a similar way as the completely optically 
thin lines. 
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For optically thick lines, a more complex dependency is seen 
in Fig. [5] (e.g., the In — Ooo line). For such self-absorbed lines, 
the wings (which are optically thin) increase with L, while any 
absorption features also become more prominent. Although the 
peak temperatures do not follow the proposed power law of the 
optically thin lines, they will become easier to detect for more 
luminous sources. 



4.2. Density 

Fig.|6]shows four lines for a range of densities, varying by more 
than an order of magnitude, corresponding to envelope masses 
ranging from 0.15 to 1.93 M . The density at 1000 AU has a 
strong influence on the water line emission. For optically thin 
lines, models with a constant p and constant luminosity show 
that a power law with an index of 0.75 accurately describes the 
dependence of line strengths on B0CH2) (see Fig. For optically 
thick lines the line will decrease in strength at higher densities, 
since the emission from the warm inner part is veiled behind a 
larger, cold water column. 

A higher density can also result in the dust becoming opti- 
cally thick. This has an effect on the higher excited lines at high 
frequencies. For example, the optically thin 331 -4o4 1893 GHz 
line (see Fig. 6, right) decreases by nearly two orders of mag- 
nitude in strength when the density increases from 0.4 x 10 6 to 
5 x 10 6 cnr 3 . 

The dust is optically thin in the low-density models and opti- 
cally thick at high density. Since the emission of this higher ex- 
citation line originates in the warm inner region, optically thick 
dust can obscure over 99% of the emission of this high frequency 
line (see § 5). 



10.0 



E 

5£ 
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10 b 
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n(H 2 ) [cm" 
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Fig. 7. The modelled integrated intensity of the optically thin 
22o-2n 1228 GHz line for Xo=10~ 4 . Intensities are shown for 
the 3 different values of Xd for a constant density of 10 6 crrT 3 
(left), and for a constant luminosity of 7 L G (right). Power laws 
with indices of 0.8 (left) and 0.75 (right) are over-plotted. In the 
right panel, at 5xl0 6 crrT 3 , this line becomes optically thick for 
Xd=10~ 6 , with its intensity dropping below the relation. 



4.3. Density profile index 

The steepness of the power law used to define the density dis- 
tribution, p, also influences water line emission. Table [6] shows 
that lines differ subtly between p-2 and p-1.5. The line intensi- 
ties of models with p-2 are lower by a factor ranging from 0.5 
to 3 than those for models with p-\.5. Absorption features also 
differ in their depths for similar models with only p varied. This 
is caused by differing temperature profiles within the Herschel 
beams. Assuming equal noCHz) at 1000 AU, the warm (T >100 
K) region is significantly smaller for a model with p-2 than for 
p=l.5. Although local densities are higher by almost an order 
of magnitude, the beam dilution is too severe to overcome such 
small regions. 

4.4. H\ s O 

The models were also calculated for the isotopologue H^O. 
Similar to the L 483 models, transitions connected with the 
ground state lines are optically thick for all models. H^O ex- 
cited lines with 7 up = 2 and 3 do have optically thick centers for 
some combinations of abundance, luminosities and/or densities, 
but are optically thin at densities below 10 6 cirT 3 . Higher exci- 
tation H^O lines are always optically thin. The best-fit power 
law slopes for the optically thin H^O lines are the same as those 
for H^O. Of course, these conclusions are limited by the optical 
depth of dust at the frequency of the line, since Hl s O lines at 
high frequencies are affected by dust. The precise frequency de- 
pends on the combination of luminosity and the density profile 
(§5). Results for H^O can be found in the online appendix. 



4.5. Line profiles 

The H2O line profiles shown in Fig. [5] and [6] all contain an op- 
tically thin and optically thick part as discussed in §3.4. Since 
the line wings of most transitions are optically thin, a 'bump' 
is present in these wings at the velocity for which that particu- 
lar line becomes optically thick. This velocity depends solely on 
the total column of water at this velocity and not on the luminos- 
ity. The height of this bump is determined by both the luminos- 
ity and density profile. The line center changes according to the 
actual optical depth of both the dust and the water itself. Line 
centers can (i) be effectively thin and show a Gaussian profile; 
(ii) be dominated by a small region and be flat-topped; (iii) show 
a self-absorption feature determined by the optical depth of the 
line only; (iv) have absorption into the dust continuum. 

4.6. Model limitations 

Our spherically symmetric envelope models have several lim- 
itations, which can affect the interpretation of observed line 
strengths. 

- The presence of a large velocity gradient, such as seen 
in infall and rotation of envelopes around high-mass stars, 
can significantly broaden the water line and simultaneously 
make it less optically thick. The velocity field in these mod- 
els was taken to be small, with an average infall velocity of 
4 km s _I at R m . For low-mass sources, such an infall veloc- 
ity combined with a turbulent width of 1 km s _1 has been 
show n to accurately r eproduce observed emission line pro- 
files dJ0rgensenl2OO4l) . 
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Fig. 6. Examples of the density dependence of line profiles for a model with p=l.5, Xq-10 , , L-l L . From left to right 

are shown the ln-0 o 1113 GHz, 2 02 -ln 988 GHz, 2 n -2 Q2 752 GHz and 3 3 i-4 04 1893 GHz lines. Densities are 4xl0 5 (solid), 10 6 
(dot) and 5xl0 6 (dash) crrT 3 . For most lines, the higher density, and thus higher water column, increases the emission, although 
absorption features are also enhanced. For the 33i-4o4 line at 1893 GHz, the dust becomes optically thick at lower temperatures (i.e. 
larger radii) within the envelope for higher densities. 



- Outflow contributions. No outflow contribution has been 
taken into account in our models. It is assumed that emis- 
sion from outflow material can be identified with sufficient 
spectral resolution in observational data with Herschel-HIFI. 
However, material heated by outflow shocks in high density 
regions will not necessarily be accelerated to high velocities. 
This material can add significantly to the emission of water 
lines, but requires a deta iled shock model to quantify (e.g., 
iKaufman & Neufeldfl996h . In addition, outflow cavities al- 
low emission from the high density regions of the inner en- 
velope to escape along paths with much lower density than 
encountered in the cold outer envelope. 

- Large extended molecular cloud or fore-ground cloud contri- 
butions can both emit rotational water lines or absorb radia- 
tion coming from the protostellar envelope. This is especially 
likely for the three transitions connected to the ground-states 
of ortho- and para-water. For example, for a cloud column 
density of 10 22 cirT 2 and an assumed water abundance of 
1CT 8 with a gas density of 10 5 cm 3 , the ground state ortho- 
H2O 557 GHz line will have an intensity of 50 mK if the 
emission fills the beam. The optical depth of this line is large, 
t =3 20, enough to absorb radiation coming from behind this 
cloud. Results from SWAS and ODIN (See § 6.2) indeed do 
not show the deep absorption in the 557 GHz line, seen in 
Fig [3|-|6] It is likely that extended cold cloud material domi- 
nates close (<0.5 km s~') to the line center in these data due 
to the much larger observing beams than Herschel. For ex- 
cited lines, emission and optical depth in a foreground cloud 
are negligible. At high (> 1 THz) frequencies the dust within 
such cloud complexes can absorb the emission from excited 
water, however. 

- A large massive gas disk can contribute to the emission of 
water. However, most gas within disks is at low (T ~30 K) 
temperatures. Most of the water will be frozen out, except in 
the inner few AU and in the warm surface layers. In addi- 
tion, high dust columns can completely obscure frequencies 
above 1 THz. Water emission from disks is thus expected 
to be orders of magnitude lower for lines at all frequencies, 
except for the transitions connected to the ground state be- 
low 1 THz. However, the small angular size of disks causes 
much of the emission to be diluted to much lower values (a 
few mK) within the Herschel beam than commonly found 
for envelopes (a few 100 mK to a few K). 



Fig. 8. FIGURE TO BE DOWNLOADED FROM WEBSITE 
Cartoon of protostellar envelope models. For two different mod- 
els, the t = 1 surfaces at 557 and 1669 GHz for the dust 
are shown. A model with low-density («o(H2)=4xl0 5 cm -3 , 
M=0.15 M ) is shown in the top panel; here the 557 GHz 
line contour is absent, since dust is optically thin throughout 
the envelope at this frequency. A model with a high density 
(n (H 2 )=5xl0 6 cirT 3 , M=1.93 M Q ) is shown in the bottom 
panel. The long-dashed line indicates the T = 100 K contour. 
In both models, the 1669 GHz line is unable to probe the inner 
warm region, regardless of its own optical depth, since the dust 
t — 1 surface lies far outside the T - 100 K radius. 



5. Effects on H 2 excitation and line formation: 
dust and micro-turbulence 

5.1. Dust 

The presence of dust in protostellar envelopes influences wa- 
ter emission lines significantly in two ways, as can be seen in 
Table [7] First, the dust can become optically thick, preventing 
water emission from deep in the envelope to escape. Second, 
the far-infrared radiation can pump the water lines. Dust be- 
comes optically thick at specific columns, depending on the 
value of the opacity, k, at a given frequency (see Fig 5c in 
lOssenkopf & Henningl (|l994) for the precise relation). Since k v 
oc v 8 , with the dust opacity index (3 typically 1.5, this is especially 
likely at higher frequencies. Table [8] column 4, shows the total 
gas and dust column densities necessary to reach a dust optical 
depth equal to 1 for each line. 

Fig. [8] illustrates this effect in a cartoon for two different 
envelope models. Two models, one with a low density (low 
M env w 0.15 Mo) and one with a high density (high M env w 1.93 
M Q ), are compared for two ground-state lines. The dust at 557 
GHz only becomes optically thick for the more massive enve- 
lope due to very high column densities, while the low-density 
model remains optically thin throughout. The t = 1 surfaces 
correspond to temperatures well over 100 K. Since the 557 GHz 
water line is optically thick at lower columns, the dust will have 
no effect on this line, except for weakening the line due to dust 
absorption. Even the H^O column is large enough to prevent 
the dust opacity from influencing the 548 GHz line. In contrast, 
the dust at 1669 GHz becomes optically thick at a much smaller 
column, which corresponds to a radius and temperature much 
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Table 7. Integrated intensity and peak brightness of lines with 
and without dust for a model with p-1.5, L = 7 Lboi and 
tt H2 = lxl0 6 cirT 3 . 



Transition 


Dust 




No Dust 






Intensity 


-^peak 


Intensity 


-^"peak 




[Kkm s- 1 ] 


[K] 


[Kkms 1 ] 


[K] 


1 io- lot 


0.95 


0.47 


1.41 


0.65 


2l2"loi 


-2.6(-2) 


0.26 


0.16 


0.39 


3l2"3o3 


0.67 


0.18 


0.39 


0.06 



further out into the envelope. The actual radius depends on the 
model in question. In the low-density model, the dust only turns 
optically thick at radii close to the 100 K water evaporation ra- 
dius, while dust in the high-density model already obscures the 
line in the outer regions of the envelope. Optically thin lines at 
this frequency will not be able to probe the warm inner region, 
since the optically thick dust effectively absorbs all water line 
emission. 

Fig. |9] shows example spectra of the 1 10-I01 line at 557 GHz, 
2i2-l i line at 1669 GHz and the 3 12 - 3 03 line at 1097 GHz for 
a model with L=7 L , n H2 = lxl0 6 cm" 3 , p=1.5, X =10 4 and 
X&= 10 4 , with and without dust. The spectra with dust are signif- 
icantly weaker by up to 50%, even for optically thin dust at 557 
GHz, since dust within the line of sight absorbs the water line 
photons. Lines that are optically thin at a frequency where opti- 
cally thick dust exists are weaker by factors up to a few or more 
if dust is present, due to less water being 'visible' . Optically thin 
lines can be pumped by far-infrared dust radiation coming from 
the region at radii larger than the r = 1 surface of the dust at the 
frequency of the line. 

For each envelope model the column density of dust depends 
on a combination of luminosity L (through the location of R m , 
Table 2) and envelope structure parameters «o and p. If the t = 1 
surface is approximated as a hard surface where all water pho- 
tons are absorbed, a typical temperature can be extracted for 
each combination of L, «o an d p. This temperature is the up- 
per limit to which water lines can probe for each model before 
dust becomes optically thick. Table[8]presents the r = 1 temper- 
atures for four astro nomically distinct mod els, based on the re- 
sults from Table 6 of |j0rgensen et alj (|2002). Model 1, with L=l 
L Q , p-1.5 and n =5xl0 6 cm 4 , corresponds to a very young 
Class YSO, which is dominated by a large, cold and massive 
protostellar envelope. A good example is the well-studied NGC 
1333 IRAS 4B source in the Perseus molecular cloud. Lines with 
a frequency above 1600 GHz are unable to probe the region 
with r>100 K. Dust is optically thin throughout the envelope 
for frequencies lower than a 1000 GHz. Model 2, with L—l L 
, p-2 and «o=5xl0 6 cm 4 , has similar parameters but its mass 
is more centrally condensed. Due to the very high column den- 
sities, dust becomes optically thick at all frequencies above 700 
GHz, making it even more difficult to probe the central warm re- 
gion than in Model 1 . Highly excited water lines are effectively 
shielded from observers within this envelope at all frequencies. 
The third model, L-2 L G , p-2, «o=4xl0 5 cm 4 , corresponds to 
a more evolved source with a low luminosity and a small enve- 
lope with most of the mass concentrated in the center due to a 
steep density profile. The relative amount of warm gas in this 
envelope is significantly higher. Such a model would correspond 
to most Class I sources, e.g., TMC-1A. The last model included 
has L=25 L , p=1.5, «o=lxl0 6 cm 4 corresponds to a higher 



luminosity source, but with a low total mass. The mass of the 
envelope is distributed over larger radii and the luminosity is 
higher. For these last two models, dust will not hide the inner 
parts of the envelope at all. Only the very warm region (T >150 
K) is veiled near the inner radius of Model 3. 



5.2. Micro-turbulence 

The structure of micro-turbulence, defined in our models as a 
velocity dispersion, during the embedded phases of star forma- 
tion is not well known and theoretical models do not yet agree 
on a correct treatment. A good rev iew of the theoretical consid- 
eration s and probl ems is given in Ward-Thompson & Bucklevl 
(1200 lL their § 3.21. |j0rgensenl(l2OO4l conclude that the inclusion 
of micro-turbulence as a velocity dispersion gives a more 
accurate reproduction of line intensities and profiles as oppose d 
to the inclusion of infall (see al so [jeirgense n et al.l d2005)). 
IWard-Thompson & Bucklevl d200U § 4.6) test different values 
of micro-turbulence, as well as different power law values 
for the variation of micro-turbulence with radius. This results 
in a large variation in the molecular line profiles of CS and 
HCO + . Similarly, Fig. [TOl shows that the choice of the velocity 
dispersion can have a large effect on all water lines, both 
optically thick and thin. All three lines (557, 1669 and 1097 
GHz) show large differences in both line profile and in total 
integrated emission up to a factor of 3, when the velocity 
dispersion is varied between 0.83 (dotted, b =0.5), 1.67 (solid, 
b =1) and 3.34 (dashed, b =2) km s _1 . Both absorption into 
the continuum and in the line are broader for larger values 
of micro-turbulence, while peak separation increases with 
increased velocity dispersion. Optically thin lines and line 
wings br oaden if b is larger. Such effects a re similar to those 
found bv lWard-T hompso n & Bucklevl d200ll) 



For molecular lines in low-mass protostellar envelopes , 
IWard-Thompson & Buckievl d200ll) and LMrgensen et all d2005h 
both adopt a constant value for the micro-turbulence with radius. 
In our case, a velocity dispersion, AV ~ 1 .67 km s 4 (b-\ km s 4 ) 
has been chosen, which best reproduces line emission for a large 
variety of molecules in sources such as L483, including CO, CS, 
but also complex molecules such as H2CO and CH3OH trac- 
ing both the warm and cold region. We recommend that for the 
modelling of water lines from specific sources, the velocity dis- 
persion for the cloud in question as derived from other molecular 
lines is adopted. 



6. Observational studies 

The conclusions presented in previous paragraphs can be com- 
pared to current observational studies. However, due to the ap- 
proximations used within the models and the unique combina- 
tions of parameters for individual sources, it is recommended 
that these models are used only as a starting point for interpreting 
water emission of protostellar envelopes from specific sources. 
As discussed in §4.6, there are several factors which will com- 
plicate any analysis of observed water lines. To compare the en- 
velope models presented above with any individual observation 
other than those from Herschel, the line intensities have to be 
corrected for the actual beam size of the telescope. 
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Fig. 9. Example spectra of the 1 10-I01 557 GHz line (left), 2i2-loi 1669 GHz line {middle) and 3i2-3o3 1097 GHz line {right) for a 
model with p-1.5, L = 7 L and n H2 = lxl0 6 cm 3 . The solid lines represent the spectra with no dust within the envelope, while 
the dotted line shows the same line with dust included in the model. The influence of dust can be clearly seen. The peaks of the 
optically thick lines are reduced in strength, the depth of the absorption changes, and the optically thin lines are significantly pumped 
by far-infrared radiation. 

Table 8. The dust properties at selected water line frequencies for four typical models. The column density 7Vh 2 required to produce 
t = 1 for the dust at each frequency is given with the corresponding radii and temperatures. 
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" Total gas and dust 

* - indicates that the dust is optically thin throughout the entire envelope 
c See text, L=l L Q , p=1.5, n =5xl0 6 cm" 3 
d See text, L=7 L G , p=2, rc =5xl0 6 cm" 3 
e See text, L=2 L Q , p=2, rc =4xl0 5 cm" 3 
f See text, L=25 L G , p=1.5, n =lxl0 6 cm" 3 




Fig. 10. Example spectra of the 1 10-I01 557 GHz line {left), 2i2-loi 1669 GHz line {middle) and 3i2-3o3 1097 GHz line {right) for 
a model with p=l.5, L = 7 L and nH, = lxl0 6 cm" 3 . The dotted line is for a model with AV = 0.83 km s _1 . The solid line for 
AV = 1.67 km s" 1 and the dashed line for AV = 3.34 km s _1 . The large effect of micro-turbulence for all lines can clearly be seen. 
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6.1. ISO-LWS 

Water has been detected in several low-mass YSOs (both Class 
and Class I) with ISO-LWS such as NGC 1333 IRAS 2/4, 
IRAS 16293-2422 Elias 29, L 1448 , L 483 and SSV 13 (e.j 
Liseau et a ll 119961: ICeccarelli et al. 



1991 iNisini et all fl99S 



Gian nini et al.1 120011: iMaret et al.l 120021) . Up to 14 water lines 



were detected, ranging from 180.5 pm to 67.3 pm. In most 
sources the total number of lines is significantly lower, but suffi- 
cient to prove the existence of water. For examp le, only a single 
line w as detected for the Class I source Elias 29 dCeccarelli et alj 
119991) . 

As an illustration , we consider the case of NGC 1333 
IRAS4A/B, for which IMaret et alj d2002l) reported the emission 
of 14 water lines together with high-J CO, OH, [O I] and [C II]. 
A model was put forward to explain the origin of these FIR lines 
using a collapsing envelo pe model without an o utflow. Based 
on the model developed bv lCeccarelli et al] dl996l) . the envelope 
and stellar parameters were derived from the best fit model to 
the 14 water lines. The water abundances in the envelope were 
found to include a jump at 100 K with a best fit of Xo=5xlO -6 
and Xd=5xl0~ 7 . A lower inner abundance is excluded by their 
model results, but a higher Xq abundance of the order of 5 x 10~ 5 
cannot be excluded. 

The physical parameters derived by lJ0rgensen et ail d2002l) 
and abundances proposed bv lMaret et alj d2002l) correspond best 
with our model with no =10 6 cm 4 , L-l L G and p = 1.5 with 
abundances 10~ 6 for Xq and 10~ 6 /10~ 7 for X^ . However, our 
predicted line s trengths differ from the observed and predicted 
line strengths in Maret et al. ( 2002), especially for low excitation 
lines. Our 22i-2i2 and 2i2-loi predictions are 0.5 10 20 and 1.1 
10" 20 W cirT 2 for the 10~ 6 /10~ 7 abundance combination within 
the ISO-LWS beam, significantly lower than the observed val- 
ues of 1 1 and 27 10~ 20 W cirT 2 . Although a inner abundance of 
10 produces brighter lines for the presented models, these lines 
are still an order o f magnitude lower than the observed fluxes in 
IMaret et all d2002l) . 

The most likely explanation is the different treatment of dust 
between the two models, which lowers the fluxes of the high- 
frequency lines in our case. However, this also suggests that our 
models can only explain the ISO-LWS fluxes if there are cav- 
ities in the envelopes through which the photons can escape. 
Recent results with the Spitzer Space Telescope show numerous 
spectrally unresolved water lines from 20 to 30 pm., which sug- 
gest that NGC 1333 IRAS 4B is indeed viewed almost face-on, 
allowing the ro-vibrational lines to es cape the protostella r sys- 
tem unhindered along the outflow axis dWatson et alj2007l) . This 
can also help the longer wavelength lines coming from the in- 
ner envelope. Another possibility is that significant contributions 
come from large-scale outflows. This is confirmed by SWAS (see 
§6.2). 

6.2. SWAS 

SWAS0 observed only the ground-state transition of 0-H2O at 
557 GHz and o-H^O at 548 GHz. Although the beam of SWAS 
(3.3'x4.5') is significantly larger than that of ISO-LWS, the in- 
creased spectral resolution of better than 1 km s spectrally re- 
solves most lines in star-forming regions. 

SWAS observed the N GC 1333 IRAS4 region dNeufeld et alj 
2000: iBergin et all 120031) where, contrary to the analysis of 
Maret et all d2002l) . they observed a wide (24.4 km s~') and 



bright ( j T* A dV = 0.5 - 1.9 K km s" 1 ) line. It was concluded that 
these lines originate within the outflow of one or m ore sources 
inside the NGC 1333 region. The model proposed by Mare t et all 
d2002l) would have produced a narrow (<4 km s~') and weak 
(0.02-0.04 K km s~') 557 GHz line for each source contained 
within the beam. Our grid gives a typical integrated intensity 
within a SWAS beam of / T UB dV = 0.5 K km s 1 for a X Q /X A 
com bination of 1 0~ 4 /10 ~ 6 with a typical line width of 6 km s _1 . 

lAshbv et al.l d2000l) report on SWAS data for B335, an iso- 
lated Jlk«s_0source with a 3cr upper limit of <0.06 K km 
s -1 . lEvans et al.l d2005l) used a wide variety of molecular emis- 
sion lines to model the chemistry and physical structure of this 
source. The low luminosity ( 4.5 Lp ) and low density (4xl0 5 
cm 4 ) found by lEvans et all d2005l) indicate that this class 
source would be represented best by a model with p=1.5, L-2 
L and « (H2)=4x 10 5 cm 3 . Our predicted 557 GHz line 
strength corrected for the much larger SWAS beam is 0.02-0.06 
K km s _1 for outer abundances from 10~ 8 to 10~ 6 , just below the 
upper limits from SWAS. The only model that would have been 
detectable with SWAS has a constant abundance of 10~ 6 . The 
high jump abundance model with Xq=\0~ a and Xd=10 -6 can- 
not be ruled out, since most of the line is self-absorbed. Similar 
conclusion holds for most other Class sources detected with 
SWAS. 

Using the SWAS archival upper limits were found on sev- 
eral Class I sources, including RCrA IRS 5 and 7. Only upper 
limits of 20-50 mK in peak temperature could be derived. These 
observed limits are higher by up to a order of magnitude than the 
predictions presented in §4 for models that would fit the physi- 
cal parameters of Class I sources. The cold envelope material is 
constrained to 10"-30", a small region within the SWAS beam. 
Thus, our models are consistent with the scarce observational 
data on Class I sources. 

In summary, both the ISO and SWAS data stress the need 
for Herschel Space Observatory data to properly constrain the 
origin of the observed H2O emission and disentangle large-scale 
and outflow contributions from the emission originating in the 
quiescent gas of the protostellar envelope. 



7. Summary and conclusions 

This paper presents several models examining the excitation of 
water and its oxygen- 18 isotopologue in envelopes around low- 
mass protostars, with parameters covering the entire range of 
luminosity and envelope characteristics of observed Class and 
I sources. A model for the Class source L 483 was used to 
illustrate the dependence of rotational water lines on the water 
abundances in the inner (T > 100 K) region Xq and outer (T < 
100 K) region X A . In turn, a larger grid of models was explored 
for the full range of parameters and likely abundances. The main 
conclusions are 

- The strengths and line shapes of all rotational water lines 
depend sensitively on the physical parameters of the envelope. 
Thus, an accurate physical model, including a well-constrained 
density profile, is a necessary pre-requisite to model water emis- 
sion and derive abundances. 

- Ground state lines from both H^O and H^O become opti- 
cally thick at low water columns. Even though water abundances 
are expected to be very low in the outer envelope (10~ 8 with re- 
spect to H2), the ground-state lines are still dominated by emis- 
sion from the cold (T < 30 K) outer (r > 1000 AU) envelope. 



See http://cfa-www.harvard.edu/swas/ for more info 



6 See http://lambda.gsfc.nasa.gov/product/swas/ for more informa- 
tion 
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Ground state line profiles are characterized by a double peak 
feature with an absorption around narrow, optically thin emis- 
sion line wings. The total integrated intensity can not be used to 
trace the water abundance. Spectrally resolving the ground state 
lines is essential if these lines are to be used as probes for the 
water abundance. 

- Higher excitation lines (E up > 150 K) are usually domi- 
nated by emission coming from the warm inner region. At very 
high excitation (E up > 200 - 250K), the lines are optically thin. 
Although these lines are weak (<0.4 K km s _1 ), they provide 
valuable information about the characteristics of the warm gas 
present in the inner regions of the protostellar envelopes. Power 
laws with an index of 0.8 and 0.75 were found to accurately rep- 
resent the dependence of optically thin lines on luminosity and 
density. 

- Hj 8 is a reliable alternative tracer of the abundances in 
both the warm (T > 100 K) and cold (T < 100 K) regions of the 
envelope. However, such lines are much weaker with predicted 
peak temperatures of only a few tens of mK, requiring very long 
integration times. 

- The presence of dust greatly influences the emission of 
most water lines. The optical depth of dust, which depends on 
the luminosity (due to varying R\ n in our models) and density 
profile parameters, can prevent optically thin water lines from 
escaping the inner envelope, resulting in orders of magnitude re- 
duction of the water fluxes. In addition, even optically thin dust 
absorbs up to 50% of the total line emission. Furthermore, the 
far-infrared continuum can pump optically thin water lines. 

- The value of the micro-turbulent velocity dispersion and 
the structure of any systematic velocity field can significantly 
influence the line profiles and intensities of all water lines. The 
adoption of a micro-turbulent width that is constant with radius, 
as inferred from line profile fits of other molecules such as CS 
and HCO + , is expected to produce the best-fitting H2O line pro- 
files, but the actual values need to be determined on a source-by- 
source basis. 

To conclude, although water is a complex molecule to model, 
total intenstities and line profiles of rotational water transitions 
can be used to constrain essential information about the proto- 
stellar envelope. 

The increased sensitivity and resolution, both spectrally and 
spatially, available with HIFI and PACS on Herschel compared 
with previous missions will provide the technology to observe 
water lines covering a wide range of excitations for the wa- 
ter molecule. The increased spectral resolution will be essential 
to locate the origin of warm gaseous water and understanding 
the structure of low-mass protostellar envelopes. Abundances in 
both the cold outer region and the warm inner region are best 
probed by a combination of lines and must include detailed line 
profiles. 
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Table 9. Integrated intensities, J TmbAV [K km s 1 ], for H2O lines in the wide parameter grid for all abundances. Ortho lines for /?= 1 .5 
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a The abundance combinations are as follows: 1: X Q = 10 4 , 2: X Q = 10 6 and a:X d = 10 6 , b: X d = 10 7 , c:X d = 10 
Note: A(B) = Axl0 B 



Table 10. Integrated intensities, J TmbAV [K km s '], for H2O lines in the wide parameter grid for all abundances, Ortho lines for p=2.0 
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a The abundance combinations are as follows: 1: X Q = 10 4 , 2: X Q = 10 6 and a:X d = 10 6 , b: X d = 10 7 , c:X d = 10 
Note: A(B) = Axl0 B 



Table 11. Integrated intensities, J TmbAV [K km s 1 ], for H2O lines in the wide parameter grid for all abundances. Para lines for p-1.5 
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6.4(-2) 


2.1(-2) 


3.4 


1.8 


0.66 3.4 


1.8 


0.60 


0.86 


0.41 


0.15 


0.87 


0.41 


0.15 
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5.0(-2) 


5.0(-2) 


5.0(-2) 


6.7(-3) 


6.6(-3) 


0.59 


0.61 


0.62 


3.6(-2) 2.6(-2) 


2.7(-2) 


0.16 


0.16 


0.16 


2.0(-2) 


1.6(-2) 


1.5(-2) 


4.3(-2) 


lxlO 6 


5.0(-2) 


5.0(-2) 


5.K-2) 


1.9(-2) 


1.4(-2) 


1.4(-2) 


0.58 


0.57 


0.57 0.14 


0.12 


0.12 


0.17 


0.15 


0.15 


5.2(-2) 


4.3(-2) 


4.3(-2) 


5xl0 6 


2.6(-2) 


1.8(-2) 


2.0(-2) 


1.2(-2) 


l.K-2) 


0.54 


0.39 


0.37 


0.43 0.27 


0.25 


0.14 


8.8(-2) 


7.9(-2) 


0.12 


6.3(-2) 


5.6(-2) 


4.3(-2) 



a The abundance combinations are as follows: 1: X a = 10 4 , 2: X Q = 10 6 and a:X d = 10 6 , b: X d = 10 7 , c:X d = 10 
Note: A(B) = AxlO B 



Table 12. Integrated intensities, J TmbAV [K km s '], for H2O lines in the wide parameter grid for all abundances for para p =2.0 



Density (cm 3 ) Integrated Intensity J Tub A V [K km s '] 
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0.34 
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0.30 
-0.26 
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0.91 
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0.44 
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0.43 
0.44 
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0.46 
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0.68 
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0.29 
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4xl0 3 
1x10 s 
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1.3 
0.66 
1.4 


0.24 
0.46 
0.99 


y.s(-2) 
0.19 
0.45 


U.40 

0.66 
1.3 


0.26 
0.45 
0.88 


9.4(-2) 
0.19 
0.47 


1.9 
3.0 
4.7 


1.1 
2.1 
5.9 


0.54 2.2 
0.89 2.8 
2.5 6.1 


1.0 
1.9 
5.2 


0.43 
0.83 
2.4 


0.79 
1.4 
1.9 


0.56 
0.94 
2.2 
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0.40 
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1.1 
0.30 
0.50 


0.19 
0.32 
0.50 


O.O^-Z ) 

0.13 
0.25 


0.34 
0.26 


0.21 
0.32 
0.43 


8.1 (-2) 
0.13 
0.26 


1.8 
2.5 
1.2 


1.3 
2.1 
4.1 


0.75 2.1 
0.96 1.9 
1.9 1.8 


1.1 
1.8 
3.5 


0.4 
0.81 
1.9 


0.58 
0.93 
0.18 


0.54 
0.78 
1.3 
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8.4(-3) 
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-3.2(-5) 


-3.9(-3) 
-l.l(-5) 
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1.6(-3) 
-2.9(-4) 


2.6(-3) 
l.l(-5) 
-2.5(-5) 


2.3(-3) 
-1.6(-4) 
-2.6(-6) 


0.46 
0.16 
-1.7(-2) 


0.45 
0.13 
-3.0(-2) 


0.42 5.7(-2) 
0.13 6.4(-2) 
-3.1(-2) 8.6(-3) 


4.0(-2) 
3.9(-2) 
-2.3(-3) 


3.8(-2) 
3.6(-2) 
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-6.7(-3) 
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3.0(-4) 
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1.5(-2) 
9.1(-3) 
7. 3 (-4) 


1.3(-2) 
6.0(-3) 
-3.9(-4) 


0.44 
0.41 
0.24 


0.40 
0.31 
5.2(-2) 


0.43 0.26 
0.29 0.30 

1.9(-2) 0.25 


0.20 
0.20 
6.0(-2) 


0.19 
0.19 

2.8(-2) 


0.12 
0.11 
3.3(-2) 


9.3(-2) 
5.9(-2) 
2.3(-3) 


9.2(-2) 
5.2(-2) 
-5.5(-3) 


9.1 (-2) 
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3.8(-2) 
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lxlO 6 
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0.14 
0.10 
4.2(-2) 


6.8(-2) 
4.7(-2) 
8.5(-3) 


2.8(-2) 
1.2(-2) 
1.1 (-3) 


0.13 
0.10 

4.16(-2) 


6.3(-2) 
4.7(-2) 
8.6(-3) 


2.4(-2) 
1.2(-2) 
1.2(-3) 


1.4 
1.8 
1.5 


0.61 
0.74 
0.55 


0.42 1.4 
0.35 1.7 
7.1(-2) 1.4 


0.60 
0.69 
0.56 


0.49 
0.34 
8.6(-2) 


0.40 
0.42 
0.32 


0.22 
0.19 
8.0(-2) 


0.12 
7.0(-2) 
7.8(-3) 


0.40 
0.43 
0.31 


0.21 
0.19 
8.3(-2) 


0.11 
6.4(-2) 
1.0(-2) 
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4xl0 5 
lxlO 6 
5xl0 6 


3.1(-2) 
1.5(-2) 
-9.9(-5) 


2.6(-2) 
8.3(-3) 
-9.4(-4) 


2.6(-2) 
7.6(-3) 
-l.l(-3) 


2.0(-2) 
1.3(-2) 
3.3(-4) 


1.5(-2) 
7.0(-3) 
-3.7(-4) 


1.4(-2) 
6.2(-3) 
-4.9(-4) 


0.50 
0.40 
6.0(-2) 


0.49 
0.34 
1.3(-2) 


0.51 0.19 
0.35 0.23 
7.5(-3) 7.7(-2) 


0.17 
0.19 

3.0(-2) 


0.16 
0.19 
2.4(-2) 


0.13 
8.2(-2) 
2.3(-3) 


0.12 
6.5(-2) 
-6.9(-3) 


0.12 
6.2(-2) 
-8.7(-3) 
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6.3(-2) 
7.8(-3) 
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4.3(-2) 
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4.1(-2) 
-2.6(-3) 
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" The abundance combinations are as follows: 1: Xq = 


10 4 , 2: 


x = io- 6 


and a:X d 


= 10 6 , b: X d = 


10" 7 , c:X d 


= 10 8 
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Note: A(B) denotes A x 10' 



